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Abstract 

We present a study of r~ — ► K~tt + tt~ij t decay using ~669 fb _1 data, collected with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. The data is recorded at a center-of-mass 
energy 10.58 GeV. The result for the branching ratio is : 

B = (3.25 ± 0.02(stat.)±%\l(sys.)) x 10" 3 . 

We also present results of the precise measurement of the branching ratio of other 3-prong decay 
modes, r _ — ► tx~ ^ tx~ v T , r~ — ► K~ K + ir~ v T , and r~ — ► K~ K + K~v T . 

PACS numbers: 13.35.Dx, 12.15.Hh, 14.60.Fg 
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INTRODUCTION 



The tau decay into 3 pseudoscalar particles has been studied since the discovery of the 
tau lepton. Moreover, decays to final states containing kaons provide a direct determination 
of the strange quark mass and Cabbibo-Kobayashi-Maskawa (CKM) matrix element |V^ S | 
[H 0, H[. Although the measurements of the r~ — > K^-k + ti^v t decay from CLEO used 
sufficiently large statistics, the overall uncertainty is as large as ~10% [4J. Furthermore, the 
branching ratio of this decay, recently measured by BABAR [3], is visibly smaller than the 
previous world average value The resonant states in r~ — > K~7r + 7r~v T decay have been 
studied only with poor statistics by CLEO 0]. The resonance study is itself of interest for 
strange spectral function measurements and leptonic CP violation studies. Theoretically, 
two intermediate resonances, the p(770) and K*°(892), are both expected to contribute to tau 
decay into R-^^v^. t~ -> K~p° (770)(-^ n + n-)v T and r~ -> A*(892)°(^ K-TT+)n-v T 
jsj. These are shown in Ufa) and (b), respectively. 

In this study, we present new measurements of the branching ratios for r~ — > tt~7t + 7i^u t , 
t~ — > K~it + it~i> T , r~ — > K~ K + ii~v T , and r~ — > K~K + K~v T decays (Unless specified, 
the charge conjugation decay is also implied throughout this paper). The measurements of 
these kaon containing 3-prong decay modes and r~ — > 7t~7t + tt~v t decay are correlated due 
to particle misidentification, so they should be treated simultaneously. 




(a) (b) 

FIG. 1: A schematic view of the e + e~ — > t + t~ in the center-of-mass system, where one r (t + ) 
decays to the signal mode (r + — > K + ir~n + i> T ) and the other r (t~) decays to the pure leptonic 
modes (r~ — > e~v T ~V e or t~ — > ii~v T u^). Each figure shows one of two possible intermediate 
resonance states in the signal mode. 

This study uses a huge amount of tau decays collected with the Belle detector at the 
KEKB asymmetric-energy e + e~ collider Qj. Since 1999, the Belle experiment has accumu- 
lated approximately 850 fb" 1 data until the summer of 2008. The Belle detector is a large- 
solid-angle magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a 
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barrel- like arrangement of time-of-flight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals located inside a super-conducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside of 
the coil is instrumented to detect K\ mesons and to identify muons (KLM). The detector 



is described in detail elsewhere [10[. In this study, we used 668.9 fb~ x data collected on the 
T(4S) resonance, 10.58 GeV, and 60 MeV below it (off-resonance), which corresponds to 
the production of 6.15 x 10 8 r-pairs. 



SELECTION OF EVENTS 

To select t~ — » A'~7r + 7r~z/ r and other 3 prong decays, r~ — > 7t~7t + tt~v t , t~ — > 
K~K + 7T~u T and t~ — > K~K + K~u T , first, tau pair events are selected. We require the 
number of tracks to be four and the sum of charges of these tracks to be zero, when the 
transverse momentum of a track in the laboratory frame is greater than 0.1 GeV/c, and the 
track extrapolates to the interaction point within ±1 cm transversely and ±5 cm along the 
beam direction. The sum of the reconstructed momenta in the center-of-mass (CM) frame 
is required to be less than 10 GeV/c, and the sum of energy deposited in the calorimeter 
is required to be less than 10 GeV. The maximum transverse momentum is required to be 
greater than 0.5 GeV/c, to reject two-photon events which include many low transverse 
momentum tracks. To reject beam-related background, we cut on the position of the recon- 
structed event vertex, requiring it to be closer to the interaction point ±0.5 cm transversely 
and ±3 cm along the beam direction. The missing mass M^ iss = (p mit — X)tr Ptr — S 7 P-y) 2 and 
the polar angle of missing momentum in the CM frame are efficient variables for rejecting 
two-photon and Bhabha backgrounds. In the definition of missing mass, p tT and p 7 are the 
four-momenta of measured tracks and photons, respectively, and pi n n is the initial CM frame 
momentum of the e + e~ beams. We require that the missing mass should be greater than 1 
GeV/c 2 and less than 7 GeV/c 2 , and the polar angle with respect to the beam direction in 
the CM frame should be greater than 30° and less than 150°. 

Particle identification is performed to select tau events containing one lepton (electron or 
muon) and 3 hadrons (pions or kaons). The magnitude of thrust is evaluated and is required 
to be greater than 0.9, to suppress two-photon and e + e~ — > qq backgrounds, where the thrust 
is defined by the maximum of (J2i l^'Pil)/ (Si \pi\) when pi is the momentum of i-th track and 
where h is the unit vector in the direction of the thrust axis - the direction maximizing the 
thrust. We require that the angle between the total momenta of the hadrons and the lepton 
momentum in the CM system should be greater than 90°, whereby the tag side lepton and 
signal side hadrons lie in opposite hemispheres, the so-called 1-3 prong configuration. Also 
the invariant mass of charged tracks and gamma clusters for each side are required to be less 
than the tau mass. Finally we require that there are no K s , 7r°, and energetic 7 on the signal 
side. Figure [2] shows the performance of the criteria used to select the 1-3 prong sample, 
where the required conditions are shown by the vertical lines in each figure. The events 
finally selected are candidates for r~ — > hih^h^Ur, where ^1,2,3 is a hadron identified as 
either pion or kaon. At this stage, the reconstruction efficiency of r~ — > K~7t + tt~v t decay 
is ~ 28%, while the dominant background modes, r~ — > 7i^7i + 7i^ir u T , r~ — > K s n~v T , 
and e + e~ — > qq are reduced to 6.0%, 1.8%, and 0.004%, respectively. The reconstruction 
efficiency of two-photon background is so small (1.6 x 10~ 4 %) that it is negligible. 

One of the most important issues in this analysis is the separation of kaons and pions. In 
the Belle experiment, dE/dx information from the CDC, hit information of the ACC, and 
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FIG. 2: Performance of event selection criteria, (a) The angle between the total momenta of the 
hadronic system on the signal side and the lepton momentum on the tag side in the CM system, 
(b) The magnitude of the thrust, (c) The invariant mass of the hadronic system. The triangular 
points show data, the dashed lines show t~ — > K~'k + ix~v t signal MC, solid lines show two-photon 
background, and dotted lines show e + e~ — ► qq background. The number of events of each sample 
is normalized to be the same. 



time-of-flight from the TOF are used to construct the likelihood for kaon (pion) hypothesis, 
L(K) (L(tt)). Figure [3] shows the likelihood ratio of kaon identification, L(K)/(L(K)+L(7r)), 
for the charged particles versus their momenta. A clean separation between kaon and pion 
is observed in the whole momentum region relevant for this analysis, where the momentum 
of kaon or pion ranges from 0.1 GeV/c to ~5 GeV/c, and the average of momentum is 
~1.3 GeV/c. As is discussed later, we choose a relatively tight particle identification (PID) 
condition for the kaon (L(K)/(L(K) + L(ir)) > 0.9) and a relatively loose condition for the 
pion (L(K) / (L(K) + L(ir)) < 0.9). With these criteria, the efficiency of kaon identification 
is ~ 73%, and the fake rate of pion mis-identification to kaon is ~ 5%. The calibration of 
the kaon and pion identification efficiencies and their fake rates are carried out by the data 
using the kaon and pion tracks in the D* + decays, D* + — > Z)°( — s- A'~7r + )7r+ sample, where 
one knows which tracks are kaons and pions from kinematics and the charges of the tracks. 
We evaluate the efficiencies and fake rates for this calibration sample and compare them to 
the Monte Carlo (MC) expectations. From this comparison, we prepare a correction table 
as a function of track momenta and polar angles, and apply them to the MC. 

This kaon identification criterion is determined by maximizing the figure-of-merit (FOM), 
when the requirements of kaon identification likelihood ratio is varied. The figure-of-merit 
is defined as: 

FOM = - = , (1) 

where S is the number of the signal (r~ — > K~tt + tt~v t ), and iV is the number of cross-feed 
background (t~ — > 7t~7t + tt~v t and t~ — > K~K + ir~b> T ). The result of FOM is shown in 
Fig. HI where one can see that the FOM is maximal with the particle identification criteria 
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FIG. 3: Likelihood ratio of kaon identification, L(K)/(L(K) + L(ir)), versus the momentum of 
the particles. The red filled circles and blue empty circles represent the real kaons and pions, 
respectively. 



used in this analysis. For the case of electron and muon identification, the probabilities to 
be an electron or muon are evaluated using the information from ECL, KLM, and other 
particle identification detectors. Also the efficiencies and the systematic uncertainties of 
lepton identification are evaluated by using the control sample 77 — > e + e~/fi + fi~. 

To reduce the feed-down background which contain K s and 7r°, we reconstruct K s and ir° 
signals explicitly. K s is reconstructed from two charged pion tracks having invariant mass 
M(tttt) within ±13.5 MeV/c 2 of the K s mass. To improve the purity of K s , the point of 
closest approach to the interaction point along the extrapolation of each track is required 
to be greater than 0.3 cm transverse to the beam direction. The azimuthal angle between 
the momentum vector and the decay vertex vector of reconstructed K s is required to be 
less than 0.1 rad. The distance between the two daughter pion tracks at their interception 
point is required to be less than 1.8 cm, and the flight length of K s transversely to the beam 
direction is required to be greater than 0.08 cm. Finally we required that the invariant mass 
of a K s candidate should not be located in the A, A, and 7 mass range, with the daughter 
tracks assumed to be electrons or pions and protons or anti-protons as appropriate. Also 
the 7r°s are reconstructed from two gamma clusters. The energy of each gamma is required 
to exceed 50 MeV for the candidates in the barrel part of the calorimeter, and 100 MeV 
for the candidates in the endcap part of the calorimeter. Also we applied a requirement on 
the invariant mass of two gammas 1/(77) to be within ±16.2 MeV/c 2 of the 7T° mass. An 
additional selection that there should be no energetic gamma whose energy is greater than 
0.3 GeV, is applied to reject remaining events containing 7r°. 
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FIG. 4: Figure-of-merit as a function of the kaon likelihood ratio value L(K) / (L(K) + L(tt)) in 
which a particle is assigned as a kaon and otherwise assigned as a pion. The uncertainties are 
evaluated from the comparison of the efficiencies and fake rates of kaon identification with that of 
control sample, D* + -> D°(-^ K~t: + )-k+ . 



EFFICIENCY ESTIMATION 



We generated ~5 million signal mode MC data with the TAUOLA based KKMC 
12| generator, to estimate signal efficiency. Also the detector response of all MC data 



sets is simulated with the GEANT [13| simulator. The TAUOLA generator uses a decay 
model for t~ — > K~7r + ir~v T , where the process is mediated by i\i(1270) and i£i(1410) 
resonances and their interference (we will call this decay model as "normal" decay model, 
in the following). To check for a possible bias of efficiency due to the specific decay model 
used, events assuming phase-space decay of r _ — > K~it + it~v r were also generated using the 
KKMC program. The evaluated efficiencies for M(KirTr) for both cases are compared in 
Fig. The relative difference of efficiencies between the normal decay and the phase space 
decay is around 1 %, so small that it is negligible. 

The average efficiencies and the fractions of the cross-feed background for all 3-prong 
decays are summarized in Table [J where we used the normal decay model for the determi- 
nation of the efficiency. One can see that the fake rate of cross-feed from t~ — > tt^it + 7t~u t 
to r~ — > A"~7T + 7r~z/ T , is very small, but since the branching ratio of r _ — > 7r~7r + 7r~z/ T is 



large, there is a substantial contamination in the reconstruction of r 
coming from misidentified r _ — > 7r~7T + 7T~v T events. 



K ix + n v T events 



BRANCHING RATIO CALCULATION 



After applying all selection criteria, the number of the events surviving as r — > 
A~7r + 7r~z/ T , t~ — > 7t~7t + 7t~v t , t~ — > K~ K + it~v T) and r~ — > K~ K + K~v T are summa- 
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FIG. 5: The efficiency of the r — > K tt +/ k u t decay as a function of M(Ktttt). Squares and 
triangular points represent the normal decay model, and the phase space decay model, respectively. 



rized in Table [TTl The possible background contaminations in these signals are from (a) 
cross-feed from the signal modes and (b) other processes such as r~ — ► 7r~7r + 7r~7r°z/ T , 
t~ — > Kg7c~u T , and e + e~ —>■ qq. The fraction of the background coming from the other 
processes is 5% to 10%, as summarized in the Table [III (3rd column) for each signal mode. 
The 4th column shows the fraction of the the main background among them: The decay 
with 7r°, t~ — ► 7T _ 7r + 7r _ 7r°z/ r is dominant one for the r~ — > 7t~7t + tt~v t and t~ — > K~7T + 7T~~i/ T 
modes, while the background from the continuum process e + e~ — > qq is dominant for the 
t~ — > K~ K +r K~v r and r~ — > K~ K + K~v T modes. 

To take into account the cross-feeds between the decay channels, the true number of 
the yield iV i truc (i = t~ — > 7r~7r + 7r~z/ T , t~ — > K~7T + 7r~v T , t~ — > K~ K + 7r~v T , and r _ — > 
K~K + K~u T ), is obtained from the following equation 

^truc = 8r j \Nf c - iVf her ) , (2) 
j 

where iVJ ec is the number of the reconstructed j-th decay mode and j\r? ther is the number of 
the background to the j-th mode coming from the other sources such as r _ — > 7r _ 7r + 7r~7r z/ r . 



TABLE I: Summary of the efficiencies and the fractions of cross-feed. 



Reconstructed 
decay mode 


Generated decay mode 


r — > 

7T~ 7T + 7T~ v T 


r — ► 


r — > 
K~K + n~v T 


r — ► 
K~K + K~v T 


r~ — > ■k~it + 'k~v t 


0.23 


0.076 


0.023 


7.3 x 10^ 3 


T~ — > K~'K Jr TX~~V T 


0.012 


0.17 


0.049 


0.023 


t~ — > K~K + ir~v T 


4.0 x 10~ 4 


4.7 x 10~ 3 


0.13 


0.060 


T~ -► K-K+K-u T 


2.8 x 10~ 6 


1.4 x 10" 4 


2.8 x 10~ 3 


0.094 
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£ij is the efficiency for detecting mode i as mode j and is the inverse of the £y matrix. 
The value of the efficiency (and migration) matrix determined by MC is given in Table 

m 

In order to determine the branching fraction, we use events of pure leptonic decays, where 



one tau decays to r — > evv and the other decays to r — > fivv 14]] (Hereafter we call such 



events { e, \x } events). The branching fraction for the decay mode % can then be written as: 

K — /vrtruc e ^V- ^r^eVu ' gr-tg /„\ 
* ~~ * ' N ' R ' ^ ' 

where N sigjefl and e eA1 are the number of {e,/x} events and the corresponding detection 
efficiency, respectively. B T ^ e -y V and B T ^-yv are the branching fractions of r — > evv and 
r — > [£vv decays, respectively, and B T ^w u = B T ^ eVu + B^^w 

With this normalization method, it is important to measure the number of { e, /x } events 
and the corresponding efficiency precisely. Figure [6] shows the comparison of the invariant 
mass of electron and muon system and the cosine of the angle between electron and muon, 
for the real {e,/x} events and the sum of MC expectation. The MC reproduces the data 
reasonably. By subtracting the backgrounds estimated by MC from the data, we measure 
the number of { e, \i } events. 



TABLE II: Number of reconstructed events (second column), the fraction of the backgrounds 
other than 3-prong cross-feeds (third column), and the main source of other backgrounds with its 
fraction to the total other backgrounds (fourth column). 





7V rec 


^yother / jy rcc 


(%) 


Main component 


in AT other 


T~ 


— ► 7T 


TT + ir~V T 


7.856 x 10 6 


10.61% 




64.3% of t~ vr 




T~ 


-» K 


~7T + TT^V T 


7.944 x 10 5 


12.15% 




34.4% Of T~ -»■ 7T 


~7r + 7r _ 7r°i/ r 


T~ 


-» K 


~K + TT~V T 


1.076 x 10 5 


6.70% 




30.4% of e+e" -» 




T~ 


-» K 


-K + K~v T 


3.162 x 10 3 


5.45% 




53.1% of e+e" -» 


99 



Figures d El El and fTUl show the comparison of the invariant mass distribution M^tttttt), 
M(K7rir), M(KKw), and M(KKK) for the data and the backgrounds estimated from the 
MC. The cross-feed background components are re-scaled using the branching ratio value 
calculated above. 

Table II III summarizes the various contributions to the systematic uncertainties on the 
branching ratio. For all cases, the uncertainty of the track finding efficiency is the most 
important source of systematic uncertainty, which is estimated from the comparison of real 
and MC data for D* — > tcD° and D° — > niiK s (K s — > tt + tt^) decay sample. The uncertainty 
of the efficiency migration matrix includes that of the kaon identification efficiency and 
lepton identification efficiency, which are important sources of the systematic errors. The 
uncertainty of the kaon identification efficiency is evaluated from the comparison of real 
and MC data for D* + — > D°7rf, and D° — > K~tt + events, while 77 — > e + e~/fi + fi~ events 
are used for the estimation of lepton identification efficienc y u ncertainty. The effect of the 
uncertainty of luminosity and cross section of e + e~ — > t + t~ [15| are rather small, because we 
used { e, /i } events for the normalization. Trigger efficiency for 3-prong modes is ~ 86 %, and 
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FIG. 6: Distributions of the {e,/i} events: (a) Invariant mass of electron and muon system, (b) 
Cosine of the angle between electron and muon. The closed circles and the solid histogram represent 
the data and the sum of the MC expectation, respectively. The open area represent the { e, fi } 
signal events, while the tau-pairs other than the { e, fi } sample and the two-photon background 
are shown by the hatched and the gray histogram, respectively. 



its fluctuation is estimated to be ~0.6 %, using the trigger simulation program for the signal 
decays. The uncertainty due to gamma veto is evaluated by using the different selection 
criteria of gamma energy. The uncertainty on the background subtraction is evaluated from 
the error propagation of the branching ratio of tau decay modes other than 3-prong modes. 
Also the uncertainty of branching ratio of the leptonic decay of tau is taken into account. 



After taking into account the backgrounds, the efficiencies and the various sources of 
systematic errors discussed above, we obtain the branching ratio for r _ — » A'~7r + 7r~z/ r 
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FIG. 7: (a) M^nir) distribution for t~ — > / k~tt +/ k~u t for data (black points). The open histogram 
is the cross-feed from r~ — > K~ix Jr ix~v T , and the dark gray histogram is the sum of all other 
backgrounds, (b) M(irinr) distribution with background subtraction. 



decay, 

B= (3.25±0m(stat.)±° o f 5 (sys.)) x 1(T 3 . 

This value agrees well with the current world average. The branching ratios for other 
3-prong decays are summarized in Table IIVI For all modes, the accuracy is improved 
significantly especially for modes including kaons, and the branching ratios are consistent 
with the previous world average values 0. But there are some discrepancies with the most 
recent precise measurement from BABAR, where the result of branching ratio of t~ — > 
K-n + n-v T decay from BABAR is B = (2.73 ± 0.02(stat.) ± 0.09(sys.)) x 10" 3 @. 
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FIG. 8: (a) M(Ktttt) distribution for r~ — ► i^~7r + 7r _ ^ r for data (black points). The light gray 
histogram is the cross-feed from r~ — ► tt~ ir + tt~ v T , the hatched histogram is the cross-feed from 
t~ — > K~ K + ir~L' T , and the dark gray histogram is the sum of all other backgrounds, (b) M(Kinr) 
distribution with background subtraction. 



CONCLUSION 

We present results from the measurement of the branching ratio for the r~ — > K~7r + 7T~v T 
decay mode, as well as those for the r~ — > 7r~7r + 7r~z/ T , r~ — > K~ K + 7r~v T , and r~ — > 
K~ K + K~u T decay modes. The results are based on the 669 fb _1 r pair sample collected 
by the Belle detector. The result for the t~ — > K^-k + -k^v t decay mode is compared with 
previous measurements in Fig. [TTJ The statistical and systematic uncertainties are com- 
patible with a recent precision measurement from BABAR, although the branching ratio of 
this study is more consistent with the previous world average values. 
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FIG. 9: (a) M{KKn) distribution for r~ — ► K~ K + tt~v t for data (black points). The light gray 
histogram is the cross-feed from r~ — > 7r _ 7r + 7r _ ^ r , the white open histogram is the cross-feed from 
r~ — > K^TT^TT^Vr, and the dark gray histogram is the sum of all other backgrounds, (b) M(KKir) 
distribution with background subtraction. 
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FIG. 10: (a) M(KKK) distribution for t~ -> K~K + K-u T for data (black points). The light gray 
histogram is the cross-feed from r~ — > tt~ tt + ir~ v T , the white open histogram is the cross-feed from 
t~ — ► if~7r + 7r - z/ r , the hatched histogram is the cross-feed from t~ — ► K~K + tt~v t , and the dark 
gray histogram is the sum of all other backgrounds, (b) M(KKK) distribution with background 
subtraction. 
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TABLE III: Summary on the systematical uncertainties. 





T — > 


r — > 


r — > 


r — ► 




TTIT'KV 


Kmrv 


KKixv 


KKKu 


Tracking efficiency error 


+3.2/-3.0 


+3.2/-3.0 


+3.2/-3.0 


+3.1/-2.9 


Efficiency matrix and PID 


+1.5/-1.5 


+1.7/-1.7 


+1.9/-1.9 


+2.3/-2.3 


Trigger efficiency 


+0.5/-0.5 


+0.5/-0.5 


+0.6/-0.6 


+0.6/-0.6 


Luminosity 


+0.1/-0.1 


+0.1/-0.1 


+0.1/-0.1 


+0.1/-0.1 


Gamma veto 


+0.8/-0.8 


+2.5/-2.5 


+1.0/-1.0 


+0.9/-0.9 


Background estimation 


+0.3/-0.3 


+2.0/-2.0 


+0.2/-0.2 


+0.3/-0.3 


Branching ratio of leptonic decay 


+0.2/-0.2 


+0.2/-0.2 


+0.2/-0.2 


+0.2/-0.2 


Total (%) 


+3.7/-3.6 


+4.9/-4.7 


+3.9/-3.7 


+4.0/-3.9 



TABLE IV: Summary of the branching ratios. 





Branching ratio 


r~ — ► ir~ir + -K~v T 

T~ — > K~7T + 1T~l' T 
T~ — > K~ K + 7T~V T 

t~ K~K + K-v T 


(8A1 ± 0.00(stat.)t° ii(sys.)) x 10" 2 
(3.25±0.02(stat.)±%H(sys.)) x 1(T 3 
(1.53 ±Om(stat.)t Q ^(sys.)) x 1(T 3 
(2.60±0.23(stai.)lo:io(sys.)) x 10" 5 



DELPHI 97 
ALEPH 98 

CLEO 99 

OPAL 00 
CLE03 03 

OPAL 04 
BABAR 08 

(PDG 06) 
(This work) 

3 4 5 6 

FIG. 11: Summary of r~ — > K~ir +/ n~v T branching ratio measurements. 



(4.90±0.80)x10 3 
(2.14±0.47)x10 3 
(3.46±0.61)x10 3 
(3.60±0.95)x10 3 
(3.84±0.40)x10 3 
(4.15±0.66)x10 3 
(2.73±0.09)x10 3 
(3.33±0.35)x10 3 
(3.25±0.16)x10 3 

x10" 3 
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